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The photonic band structures in certain two- and three-dimensional periodic networks made 
of one-dimensional waveguides are studied by using the Floquet-Bloch theorem. We find that 
photonic band gaps exist only in those structures where the fundamental loop exhibits anti-resonant 
transmission. This is also true for quasi-periodic networks in two and three dimensions, where the 
photonic band structures are calculated from the spectra of total transmission arising from a source 
inside the samples. In all the cases we have studied, it is also found that the gap positions in a 
network are dictated by the frequencies at which the anti-resonance occurs. 
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I. INTRODUCTION 

In the last decade, significant effort has been devoted 
to the study of photonic band gap (PBG) materials, 
in which the dielectric constant is periodic in spaced. 
The existence of a spectral gap in which electromagnetic 
waves cannot propagate in any direction offers the pos- 
sibility to confine and control the propagation of elec- 
tromagnetic waves. It can give rise to many interesting 
physical phenomena, as well as wide applications in sev- 
eral scientific and technical areas. For examples, on the 
scientific side, it has been suggested that the Anderson 
localization, which is complicated by the e-e interactions 
in an electronic system, is much easier to achieve inside 
the PBGs of a disordered photonic system^. In practical 
applications, ultra-low-threshold semiconductor lasers lie 
in the existence of robust microcavity modes in PBG 
systems, facilitated by the suppression of spontaneous 
emission of light 4 . Besides periodic structures, PBG has 
also been found in certain quasi-periodic systems in two 
dimensions£i£i£. 

Recently, a different kind of band gap structure for EM 
wave systems has been proposed. These are systems of 
networks made of one-dimensional waveguides*^ , which 
are experimentally realizable, even in an undergraduate 
laboratory. In particular, by using coaxial cable as the 
ID waveguide, a large photonic band gap has been found 
in a 3D network in the diamond structure, in which cer- 
tain ID waveguides have been replaced by resonant loops 
to produce strong scatteringiS. In the presence of defect 
and randomness, both defect states and the Anderson 
localized states have also been observed inside the gap. 
Furthermore, network systems are also interesting for the 
studies of many other potential physical problems such 
as PBG in fractals 11 , and the wave localization in dimen- 
sions larger than three. 

In this work, we have calculated the band structures 
of some basic lattice structures in two and three dimen- 
sions. We have also studied the photonic band structures 
in certain quasi-periodic networks in two and three di- 
mensions by calculating the spectra of total transmission 



arising from a source inside the samples. In all the cases 
we have studied, it is found that the existence of large 
gaps in a network relies mainly on the presence of tri- 
angular loops in the fundamental building blocks of the 
lattice. Each loop is capable of producing anti-resonances 
at certain frequencies with zero transmission. These fre- 
quencies coincide with the positions of the gaps found 
in the network. The outline of this paper is as follows. 
Section II describes the methods to calculate the band 
structures and the spectra of total transmission. The re- 
sults of certain periodic are also presented. Section III 
presents the results of certain quasi-periodic networks. 
Section IV gives the summary and discussions. 



II. PBGS IN PERIODIC NETWORKS 

A network consists of certain nodes connected by seg- 
ments of ID waveguides. The wavefunction at each node 
satisfies the following network equation 9 : 



-^^[cot(w%/c)] 



sin(u>lij/c) 



W = 0, (1) 



where (pi is the scalar wavefunction at node i, is the 
length of the ID waveguide between nodes i and j, uj 
and c are, respectively, the frequency and speed of the 
waves. The summation of j is over all the nodes linked 
to node i. In the case of coaxial cable, the function ipi 
represents the voltage at node i. If all the ID waveguides 
have the same length, I , it is easy to see from Eq. (I) 
the network is periodic in u> with a period liKcjl. Eq. 
(I) is equivalent to a tight- binding Hamiltonian in an 
electronic system with a zero energy eigenvalue, i.e., 



(2) 



where £j 



£)• cot kiij 



and Uj = 1/ sin kiij are respec- 
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tively, the correlated-site energy and the hopping energy. 
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FIG. 1: The unit cell for a (a) square network, (b) triangular 
network. All segments in the unit cells have the same length 
a. The lattice constants in (a) and (b) are respectively, a and 
a . 





FIG. 2: Left panel: the band structure of the square network. 
Right panel: the total transmission of the square network 
with 1369 nodes. The shaded region of the inset on the left 
panel shows the Brillouin zone. 



FIG. 3: Left panel: the band structure of the triangular net- 
work. Right panel: the total transmission of the triangular 
network with 1043 nodes. The shaded region of the inset on 
the left panel shows the Brillouin zone. 
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(3) 



Because of the discrete translational symmetry in a pe- 
riodic network, by using the Floquet-Bloch theorem, the 
wavefunctions at each node can be written as 



l C(m,n), 



(4) 



where k x and k v are the reciprocal lattice vectors or the 
Bloch wave vectors in x and y-directions, respectively. 
The indices m and n represent, respectively, the node 
indices in the x and y-directions. U (m, n) is a periodic 
function such that 



The band structures or dispersion relations lo(k) of any 
periodic network can be obtained by solving a set of cou- 
pled equations in a unit cell with the use of the Floquet- 
Bloch theorem. As an example, we consider an infinite 
network with its nodes arranged in a square lattice. Only 
nearest-neighbor nodes are connected by ID waveguide 
of length a. Fig. la shows a unit cell with four nodes 
labeled by (0,0), (1,0), (0,1) and (1,1). The period of the 
unit cell is chosen to be 2a for convenience, which is twice 
as large as a conventional unit cell for a square lattice. 
The choice of a larger unit cell leads to a band folding 
effect and will not affect the structure of photonic band 
gaps in the system. From Eq. (1), we can write down a 
set of networks equations for this unit as 



U( m +2,n+2) — t^(m+2,n) — ^(m,n+2) — t^(m,n)- (5) 

By using Eqs. (4) and (5), Eq (3) becomes 
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(6) 

where k = ui/c . For each Bloch wave vector k = [k x , /%), 
the corresponding eigenfrequencies uj(k) are determined 
by the condition that the determinant in Eq. (6) van- 
ishes, from which the band structure can be obtained. 

The band structure of the square network is shown 
in the left panel of Fig. 2, where the symmetry points 
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FIG. 4: A triangular loop with identical segment lengths a. 
Leads are connected at the input and output nodes. 



r, X and M correspond to the wave vectors (k x , n y ) — 
(0,0), (7r/2a,0), and (n/2a,ir /2a), respectively. This fig- 
ure shows that the square network does not possess a 
photonic band gap. It also shows that the band struc- 
tures are periodic in w with a period 7rc/a, which is half 
of the period shown in Eq. (1). This reason is given 
in the following. A square lattice is bipartite, i.e., it 
has two interpenetrating sublattices. For each eigenfre- 
quency of a given n = (n x ,K y ), we can always obtain 
another eigenstate with frequency u> + nc/a by changing 
the sign of all wavefunctions in one of the two sublattices. 
This can be seen easily from Eq. (6). If we make the fol- 
lowing changes ka — > ka + ir, {7(0, 0) — > — U(0, 0) and 
U(l,l) — ► — U(l, 1), Eq. (6) remains unchanged. Such a 
reduction in the period does not apply to the triangular 
network as it is not bipartite. 

To study the PBG structure in a triangular lattice, for 
convenience, we choose a unit cell that is plotted in Fig. 
lb, where one additional next-nearest-neighbor segment 
is inserted in each plaquette of a square lattice. How- 
ever, this additional segment has the same length a, as 
the nearest-neighbor ones. Thus the segments on the 
four sides of the square have to be bended and its lat- 
tice constant becomes a' = a/y2. The unit cell shown 
in Fig. lb produces a lattice that is topologically equiv- 
alent to the triangular lattice. The symmetry points T, 
X, Y and M correspond to the wave vectors (n x , n y ) — 
(0,0), (7r/2a',0), (0,7r/2a') and (n/2a' ,n/2a r ), respec- 
tively. Similar to the case of square lattice, we can use 
Eq. (4) to write down a set of equations for the disper- 
sion relations. The corresponding secular equation now 



takes the form 



-3cosA:a cosn y a' 



cos K y a' — 3cosfca coskci' cos K x a' 
-3cosfca cosKyd' 
20SK,,a' — 3cosfca 



0, (7) 



where The band structure calculated from 

Eq. (7) is shown in the left panel of Fig. 3. A large PBG 
is found with a mid-gap frequency at lu c — irc/a. Unlike 
the case of square network, the period of this band struc- 
ture is 2nc/a as dictated by Eq. (1) when all segments 
in a network have the same length a. 

The photonic band gap in a network system can also be 
studied by calculating the spectra of total transmission 
arising from a source inside the samples. This method 
is particularly useful for searching PBG when the struc- 
ture has no translational symmetry (e.g. quasi-periodic 
or fractal networks, etc.). In this method, we inject wave 
energy from a particular node inside a finite-sized net- 
work and connect all boundary nodes by leads, which 
allow waves to leak out. For propagating states, the in- 
jected waves can propagate to the boundary nodes and 
give a finite total transmission. However, if the frequency 
is inside a PBG, the injected waves cannot propagate to 
the boundary node and, therefore, give a vanishing total 
transmission. Thus, a gap in the spectrum of total trans- 
mission is a signature of a PBG. Below we use a simple 
network to show how to calculate the total transmission 
arising from a given source. 

Consider an equilateral triangular loop connected by a 
lead at each node shown in Fig. 4. The network equa- 
tions of this system have the following forms: 



ipo' — (3 cos ka)ipo + ipi + ip2 = 
ipv — (3cosfca)y>i + (p2 + (po = 
ipy - (3 cos ka)(p2 + <Po + ipi = 0, 



(8) 



where a is the length of each segment and tpy is the wave- 
function in the lead at a distance a away from the Node 
i. If the waves are injected from Node and transmitted 
through Nodes 1 and 2, the wavefunction in the lead at 
Node can be expressed as exp(ifcso) + r exp(— ifcso), 
where sq has the origin at Node 0. Similarly, the 
wavefunctions in the leads at Nodes 1 and 2 can be 
expressed as ii exp(ifcsi) and t% exp(z/cs2), respectively. 
Here r, ti and t% are respectively, the reflection and trans- 
mission amplitudes. Using these expressions, we have 
(po' = exp(— ika) + rexp(ika), (py = ti exp(ika) and 
(p2> ■ = t2 exp(ika). By substituting the above expressions 
into Eq. (8), we solve the transmission problem of the 
network, from which we find the reflection and transmis- 
sion coefficients as R = |r| 2 , Ti = |ti| 2 and T 2 — | £2 1 2 - 
The flux conservation requires that R + T± + T2 = 1. 
It should be noted that the reflection coefficient R and 
the total transmission T = T\ + T2 will not depend on 
the choice of positions 0', 1' and 2' in the leads, i.e., 
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FIG. 5: The total transmissions of the fee networks with 365 
(solid line) and 1099 nodes (dash line). 



these positions can be arbitrary. Now we apply the above 
method to a triangular network containing 1043 nodes. 
The wave is injected from a node located at the center 
of the sample. Transmission leads are connected at each 
node on the boundary. The spectrum of total transmis- 
sion is shown in the right panel of Fig. 3. A large gap is 
found in the spectrum of total transmission, which coin- 
cides with that of band structure calculation. Thus, both 
the band structure and total transmission show a large 
band gap with a mid-gap frequency at uj c — ire/ a and a 
width of Aw = 0.64o; c . The spectrum of total transmis- 
sion for the square network is shown in the right panels 
of Fig. 2. No gap is found in this structure. This is also 
consistent with the band structure calculation. It should 
be noted that there do exist dips in the spectrum, but the 
total transmission at the dip is not exponentially small 
compared to the gap in the triangular network in Fig. 
3. The dip is resulting from the small density of states 
around the corresponding frequency. 

In general, it is more difficult to find large full gaps 
in three-dimensions because its existence requires the 
pseudo-gaps in each direction to overlap with a large ex- 
tent. For the network system, PBG can still exist in 
certain networks in three dimensions. The reason will 
be given later in Section IV. Since the band structure 
of 3D networks is rather complicated, here we present 
the total transmission instead. The total transmissions 
of a fee network shown in Fig. 5 indicates that there is 
a large band gap with a mid-gap frequency oj c = -Kc/a. 
The width of this band gap is Alu — 0.8w c , which is 
25% larger than that of the triangular network. In Fig. 
5, the total transmissions of two sample sizes are plot- 
ted, where the solid line and dash line are, respectively, 
the total transmission of the fee networks with 365 nodes 
and 1099 nodes. It is seen that gap position remains un- 
changed when the number of nodes is increased 3 times. 
However, the total transmission inside the gap of the 
larger sample is lower as expected. In addition to the 
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FIG. 6: The total transmissions of the pentagonal networks 
with 1517 (solid line) and 3671 nodes (dash line). 



periodic structures presented above, we have also stud- 
ied the honeycomb network in 2D, the simple cubic, bec 
and diamond networks in 3D. However, no band gap is 
found in these structures. 



III. PBGS IN QUASI-PERIODIC NETWORKS 

It is well known that quasi-periodic structures do not 
have spatial periodicity, and thus photonic band gaps 
that appear in quasi-periodic structure may not have the 
same physical mechanism as in the case of periodic lat- 
tice by Bragg diffraction. However, several references 
have reported PBGs in quasi-periodic structures 5-6 ' 7 ' 8 . 
It was also suggested that a perfect long range peri- 
odic order may not be necessary to create a spectral gap 
in a dielectric PBG system provided that strong Mie- 
resonance scattering takes place at scattering centers. In 
the network system, we found that the existence of pho- 
tonic band gap is quite insensitive to both the period- 
icity and dimensionality of the lattice, but is strongly 
correlated with the Mie-like scattering inside a single 
loop of waveguides. To illustrate this point, we have 
studied three kinds of quasi-periodic networks. The two 
of them are two-dimensional networks and are respec- 
tively a pentagonal structure formed by Penrose tiling, 
with a five-fold symmetryiSiii, and a dodecagonal struc- 
ture formed by square-triangle tiling, with a twelve-fold 
symmetry^*. The third case is the three-dimensional net- 
work with icosahedral structure consisting of icosahedral 
symmetry that contains 15 two-fold, 10 three-fold and 6 
five-fold rotational axes. In this study, the icosahedral 
lattice is constructed by the six-dimensional hypercube 
projection methocU 5 ". The constructed structure has a 
node at the center of an icosahedron. The distance be- 
tween the nodes on the faces and the node at the center 
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5-fold Symmetric Quasi-periodic Network (3671 nodes, centre injection) 
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FIG. 7: The localized intensity distribution at frequency ui = 
0.247(27rc/a), for the pentagonal network with 3671 nodes. 

of the icosahedron is given by 

v = ^hl± a = 0.9511, 
' 2 

where a is the distance between two nearby nodes at the 
faces and r is the golden mean. For convenience, the 
nodes at the faces and the node at the center will also be 
considered as the nearest neighboring pairs, so that the 
largest coordination number of the icosahedral network 
is 12. We first consider the two-dimensional structures. 
The total transmissions of the pentagonal quasi-periodic 
networks, with 1516 (solid line) and 3671 nodes (dash 
line) are shown in Fig. 6. It is interesting to see that a 
large gap exists when sample size is 1516 nodes. How- 
ever, when the sample size is increased to 3671 nodes, 
new states appear inside the gap, which are characterized 
by the appearance of local peaks in the total transmis- 
sion. These states are found to be localized states. For 
instance, the wave intensity profile of a localized state at 
the normalized frequency Lua/2nc = 0.247 is plotted in 
Fig. 7. Due to the appearance of localized states at large 
samples, it is possible that this gap would be continu- 
ously diminished when the size of the network is further 
increased. A different behavior is found in the dodecago- 
nal quasi-periodic network. The total transmissions of 
the dodecagonal networks with 289 (solid line) and 4123 
nodes (dash line) are shown in Fig. 8. A large photonic 
band gap exists at lo c = nc/a. The width of this band 
gap is Alo = 0.6w c . This gap remains robust when the 
size of the network is increased from 289 to 4123 nodes. 
This suggests that a bona-fide large band gap can exist 
in the dodecagonal network. 

In three-dimensions, the icosaherdral quasi-periodic 
network is also found to posses a large photonic band 
gap at the mid-gap frequency lu c — nc/a, according to 
the total transmission shown in Fig. 9. The number 
of nodes of the icosahedral networks are 931 (solid line) 
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FIG. 8: The total transmissions of the dodecagonal networks 
with 289 (solid line) and 4123 nodes (dash line). 
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FIG. 9: The total transmissions of the icosaherdral quasi- 
periodic networks with 931 (solid line) and 2675 nodes (dash 
line) . 

and 2675 (dash lines). The width of the band gap is 
Au = 0.78w c , which is 30% larger than that of the do- 
decagonal quasi-periodic network, but is slightly smaller 
than that of the fee network. 



IV. DISCUSSION AND CONCLUSIONS 

It is remarkable that photonic band gaps exist in a 
number of network topologies without introducing any 
resonant loops. We summarize the intrinsic band gap 
properties of various periodic and quasi-periodic net- 
works in Table 1, where the fundamental building blocks 
and size/mid-gap ratio are listed and compared. In this 
Table, it is obvious that only the networks that contain 
triangular loops in the fundamental building block pos- 



6 



TABLE I: The summary of photonic band gaps found in various networks. uj c is the midgap frequency and Au is the width of 



the band gap. 


Networks 


Dimensions 


Periodicity 


T~l 1 j 1 1 '11" 11 1 

fundamental building blocks 


Band gap A/lj c 


Square 


2 


Periodic 


Rhombus 


Nil 


Triangular 


2 


Periodic 


Triangle 


0.64 


Honeycomb 


2 


Periodic 


Hexagon 


Nil 


Diamond 


3 


Periodic 


Hexagon 


Nil 


Simple cubic 


3 


Periodic 


Rhombus 


Nil 


bcc 


3 


Periodic 


Rhombus 


Nil 


fee 


3 


Periodic 


Triangle, rhombus 


0.8 


Pentagonal 


2 


Quasi-periodic 


Rhombus 


Nil" 


Dodecagonal 


2 


Quasi-periodic 


Triangle, Rhombus 


0.6 


Icosahedral 


3 


Quasi-periodic 


Triangle, Rhombus 


0.78 



a PBG is absence because localized states appear inside the band 
gap when the size is large. 



0.8 

a3 
g 
o 

Q- 0.6 



0.2 







/ ~- 

/ 

L / 
\ ' 


' ^-s 

>< \ 

\ 

V\ \ 
A \ \ 




i 


/ A \ 




\ ' 

V ''' 

"r" 
/ 

/ 

/ 


/ \ 1 


\ 

\ 

\ 

"\ 

\ 

\ 

\ 


i 

i 

/ 

"■■/ 

/ 

/ 

/ 


/ / \ \ \ 
// \ \ * 


\-' 
\ 

\ 

\ 

\ 










triangle 

square 
hexagon 

















0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



coa/2jtc 

FIG. 10: The transmission spectra of the triangle (solid line), 
square (dotted line) and hexagon (dash line) loops, which are 
respectively, the fundamental units in triangular, square and 
honeycomb networks. 

sess photonic band gaps. This suggests that the existence 
of the large PBGs is related to the wave interference effect 
inside triangular loops in close analogy to the Mie scat- 
tering by spherical scattering centers^. To confirm this 
observation, the transmission spectra of three loops with 
different connectivity are shown in Fig. 10. They are the 
triangular (solid line), square (dash line) and hexagonal 
loops (dotted line), which are the fundamental units in 
triangular, square and honeycomb networks. The trans- 
mission spectra are calculated by injecting waves to one 
of the leads and summing up the transmission coefficients 
at all others. From Fig. 10, it is found that although reso- 
nances appear in the square and hexagon loops, however, 
these resonances are not strong enough to produce a gap 



in any frequency in the respective networks. In addition 
to resonances, the triangular loop is capable of produc- 
ing an anti-resonance at frequency uj c = nc/a, where the 
total transmission is zero. This zero-transmission fre- 
quency is precisely the mid-gap frequency found in all the 
structures that possess gap in Table 1. The full- width- 
half- maximum of the transmission dip in triangular loop 
is 0.5uj c . This is also roughly the size of the gaps shown 
in Table. It should be pointed out that anti-resonance 
can exist only in ring-shape loops with odd number of 
nodes. The reason is as follows. For an even-node loop, 
an incident wave can propagate along two identical paths 
and interfere constructively at the node in the opposite 
side of the loop, thus, giving non-vanishing contribution 
to the total transmission. 

In summary, we have investigated the propagation of 
waves in several network structures in 2D and 3D, includ- 
ing both periodic and quasi-periodic structures. Certain 
networks have large robust photonic band gap, namely 
the triangular, fee, dodecagonal (12-fold) and icosahe- 
dral networks. The existence of such PBGs is found to 
be related to the existence of an anti-resonance inside 
triangular loops of the fundamental building units. The 
gap positions in a network are also dictated by the fre- 
quencies at which the anti-resonant transmission occurs. 
What we have studied here are intrinsic PBGs in network 
systems that do not have additional resonant loops. For 
those systems that do not have a gap, one can always cre- 
ate a gap by introducing additional resonant loops into 
the system as shown in Ref.— . The gap size as well as 
its position can also be 'tailor-made' by adjusting the 
resonant structure of the resonant loop. 
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